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O ■ Powerful relativistic jets are among the most ubiquitous and ener- 

Ph' getic observational consequences of accretion around supermassive 

O 

i-{ ' black holes in active galactic nuclei and neutron stars and stellar- 

■ mass black holes in x-ray binary (XRB) systems. But despite more 
^ than three decades of study, the structure and composition of these 

■ jets remain unknown. Here we present spatially resolved x-ray 
spectroscopy of arc second-scale x-ray jets from the XRB SS 433 
analized with the Chandra advanced charge-coupled device imaging 
spectrometer. These observations reveal evidence for a hot contin- 
uum and Doppler-shifted iron emission lines from spatially resolved 
regions. Apparently, in situ reheating of the baryonic component 
of the jets takes place in a flow that moves with relativistic bulk 
velocity even more than 100 days after launch from the binary core. 
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The jets of XRB SS 433 are well established, quasi-persistent and the best example 
of regular precession of the jet axis {IS). Extensive observing campaigns have placed 
constraints on the emission physics, geometry and scales associated with the baryonic 
component of the jets of SS 433. The generally accepted model for the thermal component 
of the jet is an adiabatic cooling model 5). This model assumes that the emitting 
matter is moving along ballistic trajectories and expanding adiabatically in a conical 
outflow, and that the temperature and matter density are only a function of the distance 
from the core. The matter is ejected (the physical processes of ejection are still uncertain) 
with a temperature of > 10® K (mean thermal energy per particle kT> 10 keV) (^- 
6), comparable to that measured for the inner regions of accretion disks. Close to the 
binary core, at distances less than ~ 10^^ cm, this jet plasma emits high-temperature 
continuum and Doppler-shifted lines in the soft x-ray band {6). The jet subsequently 
emits Doppler-shifted optical emission lines that are constrained to originate at distances 
of less than ~ 10^^ cm because of a lack of measurable offsets between blue- and red- 
shifted beams. Neither line-emitting region has been spatially resolved from the core of 
SS 433. Furthermore, the temperature of the optical line-emitting region is limited to 
less than ~ 10^ K by the absence of higher excitation hues such as He II and C III/N 
III (7). In the adiabatic model the temperature decreases with the radial distance R 
from the core as temperature TocR~^/^; thus, at distances greater than ~ 10^^ cm, the 
matter is expected to be at a temperature of < 100 K, too cool to thermally emit x-ray 
and optical radiation. Besides this thermal radiation, there is nonthermal synchrotron 
radiation, which traces the distribution of relativistic electrons and the magnetic field in 
the jet, observed in the radio band at distances between 10^^ and 10-'^''' cm {2). Further 
out, the jet is not detected again until degree-scale (i.e., at distances of ~ 10^^ to 10^° cm 
from the binary core) termination shocks are observed, deforming the surrounding W50 
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nebula (Fig. |1|, top), which is thought to be a supernova remnant associated with the 
formation of the neutron star or black hole in SS 433 {8, 9). 

Marshall et al. {6) have recently analyzed the x-ray line emission from SS 433 with 
the Chandra high energy transmission grating spectrometer (HETGS). They refined the 
adiabatic cooling model, limiting the thermal x-ray emission region to within ~ 10^^ cm 
of the jet base. They also discovered faint arc sec-scale x-ray emission on either side of 
the core, but they were not able to measure any emission lines and concluded that the jet 
gas would be too cool at such a distance to emit thermal x-rays. 

We have analysed a 9.6-ks observation from 27 June 2000 with the Chandra advanced 
charge-coupled device (CCD) imaging spectrometer (ACIS)-S. The reduced image (Fig. P 
resolves arc sec-scale x-ray jets approximately symmetrically placed on either side of a 
weaker core and aligned with both the sub-arc sec radio jets and degree-scale structures 
in W50. Analysis of the image indicates that photon pile-up (10) in the core artificially 
reduced the measured count rate there, although according to the orbital ephemeris of 
Dolan et al. (11), our observation was close to x-ray eclipse, phase ~ 0.9 (12), which may 
also have reduced the core flux. At a distance of 5 kpc, these x-ray jets correspond to 
physical scales of > 10^^ cm and soft x-ray (2 to 10 keV) luminosities of 3x10^^ to 4x10^^ 
erg s~^ each (about 3% of the average source total luminosity). 

We were able to spatially isolate and take x-ray spectra of the two jets. The spectra 
of the two lobes between ~ 0.8 and 10 keV (Fig. ^ can be fit by a bremsstrahlung 
continuum with a temperature of ~ 5 keV and Gaussian emission lines at ~ 7.3 keV for 
the approaching (east) jet and at ~ 6.4 keV for the receding (west) jet (Table 1). In 
both spectra there is an indication of additional emission, possibly also lines, at higher 
energies, but this is near to the limit of the Chandra sensitivity, and given our relatively 
low count rates, we did not attempt to fit this component. A hydrogen column density 
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of ~ 10^^ cm~^ was found for both lobes, consistent with the column density found by 
Marshall et al. (6) for the core of SS 433. A power law with a photon index of 2.1 ±0.2 fits 
the continuum equally well, and on the basis of our current data, we cannot distinguish 
between a power law and bremsstrahlung for the continuum emission. 

The significantly different line centroid energies measured for the east and west jets 
suggest that they may originate in the moving jet material. To test this hypothesis, 
we developed a model that sums the emission from ten equally spaced phase bins (13) 
around the precession cycle, calculating for each bin the Doppler shift 6 (14) and Doppler 
boosting of the line peak intensity 5^ (assuming that the jets are intrinsically symmetric 
about the precession axis and east to west). A simulation based on this model for a 10- ks 
observation, assuming a rest energy for the line of 7.06 keV, corresponding to Fe XXV 
k/5 and using the same normalizations for the Gaussians as found in our observation, 
corresponds well to our data (Fig. ^ and Table 1). 

The line energies and widths we obtain in our simulation are a good match to those 
measured in our observation (Fig. ^ and ^). Thus, we infer that the jet is emitting around 
the precession cone at a large distance from the core. Our model predicts twin-peaked 
lines, due to the distribution of Doppler shifts as a function of time, with a stronger 
blue peak (due to Doppler boosting); however, this twin-peak is not seen in the 10-ks 
simulation, which have a relatively low signal-to-noise ratio. However, if our model is 
correct, then additional observations should reveal this structure. Although we cannot 
use our model to conclusively establish that the Doppler-shifted line is Fe XXV k(3, we 
do need a transition of hot, highly ionized Fe in order to achieve a high enough energy for 
the line in the east jet. We cannot exclude an origin for the ~ 6.4 keV line in the west 
jet from neutral Fe; however, we cannot envisage a plausible geometry whereby we are 
observing neutral Fe at 6.4 keV in the west jet and the same line extremely blue-shifted 
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in the east jet. A further prediction of this model is that the east jet should be stronger 
than the west jet as a result of the overall larger Dopplcr boosting. The ratio of 2 to 10 
keV fluxes (line and continuum) from the east to west jets is ~ 1.3, approximately what 
would be expected for the model considered here. Although we cannot explicitly show 
that the source of the continuum emission arises in the moving jet gas, the measured 
bremsstrahlung temperature is consistent with an origin in the same plasma as the high- 
excitation Fe hues. We therefore conclude that we observed lines, and probably also 
continuum, from a hot (> 10''K) plasma which is moving with a similarly high velocity 
to that measured for the "inner" jets {6), ~ 0.26c (where c ~ 3 x 10^° cm s~^ is the speed 
of hght), on scales of 10^^ cm from the binary core. 

This is contrary to the adiabatic cooling model {4^6) which at that distance predicts 
that the baryonic component will have cooled to the ambient interstellar medium tem- 
perature; therefore, the jet has to re-heat itself. The most plausible energy source is the 
kinetic energy (~ 10^^ ergs/s) associated with the bulk motion of the jet (6). Because we 
observed optical emission lines coming from matter with temperatures of < 10^ K at a 
distance of ~ 10^^ cm, this re-heating must occur between ~ 10^^ and 10^^ cm from the 
core. 

Alternatives to this scenario can be considered. A power law flts the continuum well, 
and by analogy with the extended x-ray jets seen in active galactic nuclei (AGN) (i5), we 
cannot exclude an inverse-Compton or synchrotron origin for the continuum component. 
However, because both of these emission mechanisms require a population of hot electrons, 
they also imply reheating at those large distances from the central source. Regardless of 
which form of reheating is associated with the continuum spectrum, the emission lines 
require a baryonic component. Reflection of the core emission does not seem possible 
because the line flux we observed at 10^^ cm is only a factor of ~ 10 less than that 
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observed to originate at 10^^ cm (6), and yet our images indicate it does not subtend 
a significantly larger solid angle than the jets on smaller scales. Furthermore, the line 
to continuum flux ratio we observed is an order of magnitude greater than in the core 
[-^lineZ-^continuum is ~ 1 for our cast jct, but is ~ 0.1 for the 7.3-keV line in Chandra HETGS 
observations (6)]. The lines from these extended jets are still an order of magnitude weaker 
than the moving core x-ray lines (5, 6), and would not be expected to show measurable 
changes in their Doppler shifts (unless, e.g., a major outburst from SS 433 resulted in 
enhanced emission at some particular precession phase). Nevertheless, they should be 
present as weak features in spectra that do not spatially resolve the core and jets. 

Reheating (in situ acceleration) of a leptonic (electron and possibly also positron) 
component of jet plasma, presumably by conversion of the jet kinetic energy, is commonly 
invoked to explain synchrotron or inverse Compton energy emission on large scales in the 
jets of AGN (15-17). The reheating of the baryonic plasma that we observed in SS 433 
requires a similar tapping of the bulk kinetic energy of the flow, by processes that are able 
to act on atomic nuclei. In the internal shock model for blazars [16), a slightly varying 
jet speed results in shocks downstream in the jet flow. In the context of this model, 
random discrepancies of about 15% in the velocity of the jets from SS 433, inferred from 
long-term optical monitoring [3), could produce shocks a few hundred days downstream 
(for a precession period of 162 days and a mean jet velocity of 0.26c), comparable to what 
we observed (18) . 

To conclude, these observations spatially resolve a line-emitting region in a relativistic 
jet. They reveal that the hot x-ray emitting region is still moving relativistically on phys- 
ical scales that are orders of magnitude larger than previously inferred, with no evidence 
for substantial deceleration, supporting models for extended x-ray jets in AGN that re- 
quire bulk relativistic motion on large physical scales (15). Furthermore, this observation 
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demonstrates that particle reacceleration in a relativistic jet, previously inferred only to 
act on the leptonic component, can act also on atomic nuclei. Comparison of this highly 
super-Eddington, mass-loaded jet with those of other x-ray binary systems and AGN 
can provide unique insights into the matter and energy content of relativistic jets and 
hence the coupling of accretion and outflow in conditions of extreme gravity, pressure, 
and energy density. 
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Table 1: Comparison of observed line centroid energies E and widths a with simulations 
based on symmetric jets integrated over an entire precession cycle and a rest wavelength 
of 7.06 keV (Fe XXV kp). 
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Figure 1: (Bottom) Image of SS 433 with Chandra ACIS-S, rebinned to one-quarter 
native pixel size and smoothed with a 0.5-arc sec Gaussian (comparable to the resolution 
of the telescope). The predicted projection on the sky of the jet precession cycle is 
indicated. The blue line refers to the jet that, for the most of the time, is approaching 
Earth (it is receding only ~ 16 days over the 162-day precession cycle), and the red 
line refers to the (mostly) receding jet. The x-ray emission lies in the jet path and, 
while concentrated in two lobes, shows weaker emission extended beyond these to larger 
distances at the predicted position angles. (Top) Comparison with the surrounding W50 
radio nebula (9) shows the correspondence in alignment of the structures, both of which 
reveal the signature of the jet precession (19). DEC, declination; RA, right ascension. 
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Figure 2: (Top) Raw (i.e., not regridded) Chandra ACIS-S image of SS 433 with the 
projected precession cycle of the jets superimposed on it, indicating the regions of the 
jet we isolated and analyzed. (Bottom) The two spectra of the lobes (crosses), east 
(left) and west (right), with the fitted model [absorbed bremsstrahlung (dashed line) and 
Gaussian emission line (dashed-dotted line)] and residuals (x)- The error bars represent 
the la confidence interval for each measurement; the residuals are plotted in units of a. 



12 



P = 162.4 days 
v=0.26c 




10 1 

channel energy (keV) 



10 



Figure 3: Simulated spectra (crosses), integrated over one entire precession cycle and fit 
for the same integration time, with the same binning and model [absorbed bremsstrahlung 
(dashed line) and Gaussian emission line (dashed-dotted line)] as for our observation in 
Fig. 2. Residuals (x) of the fit are also plotted. The error bars represent the Icr confidence 
interval for each measurement; the residuals are plotted in units of a. Parameters for the 
jet precession model (jet velocity v = 0.26c, opening angle 9 = 21°, precession axis angle 
i = 78°, period P = 162.4 days) are from {3). For each of ten precession phases, the 
observed line energy E^hs = ^E^est and Doppler-boosted peak intensity Jobs = ^'^Irest were 
calculated and summed to produce the resultant line profiles (the model can be visualized 
as simulating the integrated line emission of the shaded region in the schematic). The 
normalization of the summed lines was chosen to be the same as that in our observations, 
and the rest energy of the line was chosen to correspond to that of Fe XXV kj3 (i.e., 7.06 
keV). 
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